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ABSTRACT 

An algorithm for the solution of sequence-dependent routing 
problems is presented and programmed in FORTRAN IV for use on digital 
computers. Solutions, computation times and iteration requirements 
are summarized and discussed for eleven test cases. 

With specific modification of the input data, a typical traveling 
salesman closed-loop problem may be solved by the same program. 
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LIST OF SYMBOLS AND ABBREVIATIONS 



X - a subset of all feasible solution vectors 

Y - a subset of X 

Y - the complement of Y with respect to X 

W(X) - a bound on the objective function for all possible solution 
vectors in X 



leg k - one of the sequence of arcs which form a complete route 

(the k-th leg of a route between N nodes is that arc (i,j) which 
is traversed between the k-th and (k + l)-st nodes visited 
in sequence on the route) 



arc(i,j) - a directed path from node i to node j 
]<^ 

A / (a..) - the matrix of costs of traversing arc (i,j) on the k-th leg 

of the route 

^k' ~ current working matrix of costs of traversing arc (i j) 

on the k-th leg of route. (Initially = A^ but is 

changed by the operations of the algorithm) 

a_ summed over the set of (ij:k) for committed arcs and legs 

M' - the reduced form of M, 
k k 

q(ik/3k‘*^) “ the reducing constant for 



9 (i J :k) 
P P 



the second smallest element in M,' 



0 (i ,\ :k ) = max 0 (i A :k) where k is uncommitted 
o o o , p p 



X - represents plus infinity as a matrix element 
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I. INTRODUCTION 



The algorithm programmed in this thesis, presented by DeHaemer 

[Ref. 1], uses the branch-and-bound technique to find the optimal 

route between N nodes. It determines the beginning and ending nodes 

and passes through each node exactly once. The criterion for optimality 

is to minimize total cost in traversing the (N-1) arcs of the route where 

the cost of traversing each arc is a.., which is a function of the k-th 

ij 

position in the sequence of arcs forming the route. 

The purpose of this paper was to construct a computer program 
which would solve the general class of sequence-dependent routing 
problems using the above mentioned algorithm, given the matrices of 
all possible costs for each leg of the route. The difficulty in solving 
this class of problems has been in finding a method of selection of 
tours which avoids evaluation of all the (N-1) 1 possible tour costs in 
determining an optimal route. 

Although several algorithms for typical traveling salesman 
problems have been proposed and programmed for a computer [Ref. 2], 
this paper presents the first program and results using the algorithm 
presented in the next section. 

The operational results of solving several test problems are 
given along with a discussion of the limitations of the computer program. 
It is assumed that the reader is familiar with the branch-and-bound 
technique. References 1 and 3 discuss general background of branch- 
and-bound methods . 
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TYPICAL SEGMENT OF TREE 



W(0O = 7 




Notation : 

W(X) - a lower bound on objective function for all possible solution vectors, 
attached to base node X of tree 

Y - right-hand nodes with notation as follows (i,j:k) (i.e., k-th leg 

of route from i to j) 

W(Y) - lower bound associated with node Y 

Y - left-hand nodes with notation as follows (i, j:k) (i.e. , k-th leg 

of route is not from i to j) 

W(Y) - lower bound associated with node Y 

Note: For computer application, right-hand nodes are labeled with 

positive numbers and left-hand nodes are labeled with negative 
numbers . 



Figure 1 



n. THE ALGORITHM 



The basic method employed by the algorithm is the branch-and- 
bound technique. The set of all possible routes through N nodes is 
broken up into smaller and smaller subsets and a lower bound on the 
cost of the best route in the subset is obtained. The bounds are then 
used as guides in determining further partitions into smaller subsets 
until the algorithm eventually isolates one or more subsets which are 
complete routes whose costs are less than or equal to the lower bounds 
for all other subsets. These routes are then declared optimal. 

The algorithm generates a tree whose nodes represent subsets of 
routes as illustrated in Figure 1 . The base node of the tree establishes 
an absolute lower bound on all possible routes. Each branch or segment 
of a branch is a complete route or subset of a complete route respectively. 
An example tree for an entire problem as generated by the computer 
program may be seen in Appendix A. 

It is assumed that the set of matrices A^^ can be specified for all 
(N-1) legs of the route. A problem with N nodes requires that (N-1) 
legs of a route be determined. Each leg k of a route is specified as 
being an arc (i,j) which is a directed path from node i ^ node j. 

The algorithm as used for the computer program is listed here in 
complete detail. The first three test cases in Section V. A. are worked 
out in some detail in Ref. 1 and sufficient background of the algorithm 
may also be found in the same reference. The only modifications made 
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here in this algorithm are in the branching rule of step 8, elaboration 
of step 7 for the computer program ^ and in the branching to step 7 from 
step 4 when sufficient legs of a route are known so that a complete 
route may be specified. 

The Steps of the Algorithm 
Step 1: 

The initial setup of the algorithm is made as follows: 



1 . 
2 . 
3. 



Set for k = 1,2, . . . , (N-1). 



X is the set of all possible routes. 

Set Z = oo and Leg = 0. Z will be the cost of the 
o o 

optimal route at the end of the algorithm. 



Step 2: 

Find the minimum element in each matrix and reduce the 
matrices. An absolute lower bound on the cost of all 
tours is found. 

1 . For each leg k, k = 1,2,..., (N-1) , find i . j . and 

k k 

]i^ 

q(i /j :k) such that q(i ,j :k) = min min m... 

K K. K K • • 1 J 

1 J 

‘ 'k k 

2. Reduce M, to M, where m.. = m . - q(i, / j, :k) for 

k k ij ij k k 

all i, j / and k . 

3. Label node X with W(X) = ^q(ij^/j^:k) summed over 
k = 1,2..., (N“l) . This label is the absolute lower 
bound on the cost of all tours. 
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Step 3: 



Choose the subset for the next tree extension as follows: 

'k 

1. 0 (i J :k) = min m.. for each k where leg k is 

P P . . V • • 

uncommitted . 



9 (i /j :k )= max 9 (i ,j ;k) where k ranges over the 
o o o ^ P P 



uncommitted legs . 



3. Then Y = (i ,] :k )and Y = (i :k ) are the next 
0-00 o o o 

branches fromX. 



Step 4: 

Label Y by W(Y) = W(X) + 9 (i J :k ) . 

o o o 

Step 5:^ 

Since an arc is to be committed to a leg, a new set of 
restricted matrices are formed by the following actions: 

1 . Delete M, . 

k 

o 

2. a. Delete all elements in M, , except row j . 

k + 1 o 

o 

I 

b. Delete columns i and j in M, , , . 

o o k + 1 



Step 5 of the algorithm was accomplished in the computer 
program through the use of the variable matrix ARCCOM and the variable 
DEL which allowed only certain matrices and certain elements in these 
matrices to be considered in the succeeding steps. 
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3. a. Delete all elements in M, , except column i . 

k - i o 

o 



b. Delete rows i and j in M, , . 

o o k -1 
o 



4. Delete rows i and i and columns i and j in all M, 
o o o "o k 



except in M, , , and M, 

k + 1 k - 1 

o o 



5. Relabel the matrices as M, . 

k 

6. Leg k is now committed to arc (i ,] ) . 

o o 



7. If (N-3) legs have been committed, go to Step 7 . 



Step 6: 

Initiate procedures to determine what the next leg of the 
route should be. 

1 . For each k where leg k has not been committed to a 

route, find i^, j^, and q(i^,j^:k) such that q(i^,j^:k) = 
k 

min min m. . . 



2 . Reduce to for those legs k which are not 



committed and for all i,j of uncommitted arcs where 



•k 
m. . 
ij 



m.. 



3 . Label Y by W(Y) = W(X) - summed over k 



for uncommitted legs. 
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Step !■} 

Ascertain whether a route has been determined and if it 

has an upper bound which is equal to or less than Z . 

o 

1. Increment leg by one since a leg has been committed. 

2. If (N-2) legs of route have been committed and W(Y)^ 

Z , go to Step 10 . 
o 

3 . If (N-2) legs of route have been committed and W(Y) > 
Z^ , go to Step 8. 

4 . If (N-2) legs of route have not been committed and 
W(Y) ^ Z^, go to Step 8 , substep 4. 

5. If (N-2) legs of route have not been committed and 

W(Y) > Z , go to Step 8 . 
o 

Step 8: 

Determine the node X from which to branch as follows: 

1. Make the last Y node non-terminal since it is either 

the end of a complete route or the end of a segment of 

a complete route which has a cost which is greater 

than Z . Therefore, a search of Y nodes for suitable 
o 

branch points must be made. Go to substep 2. 



2 

Note that when (N-2) legs of route have been committed, the 
last leg is automatically determined and hence computation ends when 
(N-2) legs are known. 
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2. Choose the lowest numbered left-branch node with a 

label W(Y) ^ Z and branch from this node X. Go to 
o 

Step 9. For all Y nodes with labels W{Y) > Z , 

o 

consider them non-terminal since they would all lead 
to higher cost routes. If there is no Y node which is 
a candidate for branching, go to substep 3. 

3. All nodes have been made non-terminal by substeps 1 
and 2 of this step and hence the optimal route has been 
found. STOP . 

4. If substep three of Step 7 was satisfied, make last Y 
node to be the node X from which to branch. Make Y 
node non-terminal and set W(X) = W(Y) . Go to Step 3 . 



Step 9: 



Set up the cost matrices and label node X as follows; 

1. Set leg = 0. Then determine number of legs committed 
on limb of tree from which branch is to occur and set 
leg = to the number of Y nodes on the limb. 

2. Compute g = ^a^. summed over the set of (i,j:k) for 
committed arcs and legs at this point in the tree. 

3. If no legs have been committed, set otherwise 

I 



ser M, “ • 



4 . Carry out substeps 1 thru 4 of Step 5 for each of the 
committed arcs and legs . 
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5. Block paths which are not allowed (i.e., those left- 
hand nodes encountered on this branch of tree are 
forbidden nodes) . 

6. Carry out Step 6 substeps 1 and 2 . 

7. Label X with W(X) = g summed over k for 

the uncommitted legs. 

8 . Go to Step 3 . 

Step 10: 

Determine complete route which has been found. 

1 . Arrange the committed arcs and legs to determine 
missing leg and arc on this leg. 

2 . Make last Y node non-terminal since a route was 

determined . 

3 . Set Z = W(Y) . Go to Step 8, substep 2 . 

o 

End of Algorithm 

A flow chart of the algorithm is in Figure 2 . 
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FLOW CHART OF ALGORITHM 



^ START ^ 



each matrix and reduce M)c 
to 

W(X) = sum of reducing 

constants 

Label X with W(X) 



Step 3 



I 



Choose subset for next tree 
extension by finding 

Q(ip/ 3p:h) 

Y—(io/jo'^o) (io*3o‘^o) 



Step 4 



I 



La^l Y by 

W(Y) = W(X) + 9(ioJo:ko) 



Step 5 



Form new set of by 
restricting arcs and legs 
previously committed in 
subset under consideration 



Step 6 



Find minimum element in 
new set of Mk and reduce 
Mk to Mk 
Label Y by 

W(Y) = W(X) + sum of 
reducing constants 



Step 8 



step 1 > 




Set Ak = Mk for all k. 

X is the set of all routes 


Zo = 00 Leg 


= 0 


Step 2 >j 




Find minimum element in 




Step 9 \ 


L YES 


— * 


Form the set of reduced] 
matrices mJ. for the 
branch under consider- 
ation 




Step 10 






Set Zo = W(Y) Complete 
route and obtain route 

cost. Make note of 
best route . 








YES 






Figure 2 



16 



W H O CL, 



'i 



III. PROGRAMMING CONSIDERATIONS 



The first decision that had to be made before programming of the 
algorithm began was what computer language would be most appropriate. 
Since one of the primary purposes of this project was to explore the 
feasibility of computerized solutions using the algorithm rather than to 

3 

develop an efficient program for large-scale problems, FORTRAN IV was 
chosen as the language due to its ease of application. 

One of the important factors to consider for computer applications 
is requirement for storage space. The strategy used for selection of 
the branch point in Step 8 can have a direct effect on storage require- 
ments . There are two basic strategies which may be used: 

Strategy 1 : Branch from the lowest bound. This strategy is the 
one used in the original algorithm [Ref. 1] and has the advantage that 
the total computation required to reach optimality is minimized in the 
sense that any branching performed is also that which must be performed 
under any alternate policy. Its primary disadvantage is that no terminal 
nodes are discarded and hence storage requirements may become ex- 
cessive. In addition, it brings Step 9 of the algorithm into play more 
often which requires time to backtrack through the tree and set up the 
matrices for a further branch from the chosen node. 



3 



Large-scale here is considered to be when the number of nodes, 



N, is greater than 20. 
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Strategy 2 : Branch always from the latest Y node if a complete 



route has not been determined and discard nodes from storage that are 
no longer in contention for branch points or for the optimal route. This 
is known as a "branch to the right" policy. It has as its primary 
advantage that the amount of computer storage required is minimized 
since nodes are discarded when they are no longer required. Also, 

Step 9 of the algorithm will not be called upon as frequently as under 
Strategy 1 . 

Strategy 1 was originally employed, but for the few test cases 
considered, the number of iterations and time required to obtain the 
optimal route was in general greater than that required under Strategy 2 
and hence the program presented uses Strategy 2 . 

As mentioned in Reference 1, a very useful feature of this routing 
algorithm is that one can stop at any point after the first complete route 
has been determined and have a feasible tour, although it may not be 
optimal. In the computer application of the algorithm, it may be the 
case that sufficient storage space or time required to reach the optimal 
solution may not be available. Hence, if one is willing to accept a 
suboptimal solution such as a solution below a given cost, this given 
cost could be input to the program and as soon as a solution that has a 
cost less than this amount has been found, computation can be halted. 
This may be found to be extremely useful when dealing with large-scale 
problems where to pay for sufficient computer time to reach the optimal 
solution might be prohibitive [Refs . 1 and 5] . Note that in test problem 
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Number 11, a solution within 4% of the known optimal solution was 
obtained in a very short period of time, but that nearly 300 minutes 
and 25,000 iterations later, the same solution was found and the optimal 
route had still not been located . 

Although Reference 1 gives a modification to the basic algorithm 
for symmetric matrices, the modification was not incorporated in the 
program presented here. 

For test problems 1-10 presented in Section V, storage require- 
ments did not become excessive as will be discussed in more detail 
under Section VI on computational results. 
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IV. THE COMPUTER PROGRAM 



The computer program is entirely integer in nature except for the 
variables used in conjunction with the timing routine. A detailed 
description of the major variables used in the program may be found in 
Appendix D. Originally, the program was compiled using the FORTRAN 
G“level compiler and consisted of a main program where the major 
portion of all iterations was accomplished, and two subroutines, one 
used for Step 5 of the algorithm and the other for output. 

It was noted that the FORTRAN H -level compiler generated an 
object code which was superior to the G-level compiler, particularly 
for extensive looping and arithmetic operations which were present in 
this type of program. An attempt was made to compile the identical 
program using this H -level compiler but the program size combined with 
its complexity was too large for the compiler to accommodate. At this 
point, the program was broken up into a main program and eight sub- 
routines, all of which the H -level compiler could handle. The computer 
program flow along with a brief description of the subroutines is illus- 
trated in Figure 3 . 

Maximum storage utilization was attained by specifying that 
nearly all variables be INTEGER*2 . This meant that the principal 
iteration information for the tree which was maintained for purposes of 
being able to branch from any node was limited to numbers less than 
or equal to 32,767. This limitation applies to bounds on nodes and 
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number of iterations; hence node numbers, since node numbers are 
directly related to iterations. All right-hand nodes are labeled by 
positive numbers which identify them with the iteration on which they 
were obtained and likewise, left-hand nodes are labeled by negative 
numbers . 

Another limitation of the program as presented is that the number 
of nodes be equal to or less than 2 0. The number of tours which can 
be expected to be obtained is limited to 30. The first tour is obtained 
by branching to the right immediately until a complete tour is specified 
which takes place on the (N-2)-nd iteration. All future tours must have 
cost equal to or less than the previous tour or they are not considered 
or counted as a tour for the purposes of the program. 

All of the limitations discussed are limits of the computer program 
as presented and may be easily modified by changing the appropriate 
DIMENSION statements. Iteration information contained in the matrices 
YTAB and YBTAB which is used for constructing the branch point becomes 
the primary storage-limiting factor when the number of iterations is 
expected to be in the thousands. For 150 iterations, which was used 
for the first eight test problems, the entire program required 114,000 
(114K) bytes of storage. Each increment of 100 iterations above the 
150 used requires 1 . 8K bytes of storage and therefore 2500 iterations 
as used for test problem Number 9 required 42K more bytes which led 
to a program size of 156K. 
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For the typical traveling salesman problems discussed in the 
next section, the optimal route as expressed by the computer output 
has been adjusted to reflect the actual route which excludes the dummy 
node (N+1) . Typical computer solution output for both a sequence- 
dependent case and a typical traveling salesman case may be found 
in Appendix B . 

A timing routine used in Reference 4 is included in the program 
for purposes of obtaining actual problem solution times which excludes 
all input and output buffering times. 

The program follows the algorithm step by step. Documentation 
is interspersed throughout to enable a casual reader to understand the 
basic program flow. The entire program may be found in Appendix E. 
Appendix C contains the make up of the computer card deck. 

In order to provide dynamic allocation of storage space based 
upon the number of nodes in a given problem and the number of iterations 
desired, modifications to the basic program presented in Appendix E 
have been provided in Appendix F. Details on the specific changes are 
given in Appendix F. The primary advantage of these modifications is 
that the user does not have to change all of the variable specifications 
and dimension information cards in the 8 primary routines each time 
different values for N and ITS are used (N is the number of nodes; ITS 
is the number of iterations desired) . Only the appropriate job control 
language (JCL) card which specifies the storage and time requirements 

for the execution of the program must be changed. 
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V. TEST PROBLEMS 



A. SEQUENCE DEPENDENT PROBLEMS 

The first three test cases were problems whose description places 
them into the class of sequence-dependent routing problems. Problems 
1 through 3 were taken directly from DeHaemer [1]. Problems 1 and 2 
consisted of matrices which were asymmetric. In problem 3, all 
matrices were symmetric. As was noted in Section III, the computer 
program does not provide for special treatment of symmetric matrices, 
but it was desirable to include symmetric matrices as test problems. 

Problem No . 1 



Suppose an itinerant salesman must be routed so that his 
travel expenses are minimized v/hile visiting 5 different cities. He 
must complete a leg of his route on each of 4 consecutive days. Travel 
expenses vary as a function of the day on which the travel occurs. At 
certain times, no public transportation is available and the costs reflect 
the price of the available charter transportation. All possible costs 
have been tabulated for each of the 4 traveling days and are presented 
in Figure 4 . 
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Problem No. 1: Initial Set of Cost Matrices 



Figure 4 
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Problem No . 2 



This problem has the same framework as problem 1 except 
that there are 6 different cities and thus there are 5 legs of the route. 



The matrices of all possible costs are tabulated in Figure 5. 
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Problem No. 2: Initial Set of Cost Matrices 



Figure 5 



Problem No. 3 



Figure 6 contains a set of four symmetric cost matrices 
from which a minimal cost route is desired. 
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Problem No. 3: Initial Set of Cost Matrices 



Figure 6 



These first three examples are discussed in more detail along 
with sample calculations in Ref. 1. 

It would have been desirable to have larger test problems for 
which an optimal route was known. In order to avoid the lengthy hand 
computations involved in the setup and solution of a larger problem, 
it was thought that the typical closed-loop traveling salesman problems 
which have known optimal solutions and are abundant in the literature 
could provide additional test cases. 



B. TRAVELING SALESMAN PROBLEMS 

By appropriate modification of the input data, the typical closed- 
loop traveling salesman problem (hereafter referred to as TSP) can be 
solved by the program. It was necessary that the problem be structured 
in a manner such that the route would be closed as opposed to the open- 
ended route determined by the algorithm, visiting each node exactly 
once. Since the optimal route in a TSP is independent of the starting 
node, it was observed that the addition of one dummy node and hence 
one dummy leg attained the desired results. This can best be illustrated 
by an example . 
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Suppose the following matrix of costs between 4 nodes was 



given: 

12 3 4 

1 X 2 6_ 3 N = 4 

2 5 fx”l 4] 

3 2 • 3 X 5l No . of legs = N - 1 = 3 

4 ll_7__2_x| 



It is assumed, as is usually the case in the TSP, that the matrix is 
the same for each leg. Consider the following set of four matrices 
which have one additional dummy node (node 5) besides the original 
4 from above. 



Mj M3 
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1 


2 3 
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2 3 4 5 
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X X 
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X X X X 
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X 
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fx f 


T1 


X 
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ix 1 4^ X 
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X 
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X 
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1 3 X 


5 I 


X 
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!3 X 5 • X 
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2 
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|_7__2_ 
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[7_2 xji X 
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No. 


of legs = 
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X 
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X X 


X 


X 


X 


X X X X 


X 


X 


X 


X 


X 


N-1 


= 4 



The number of nodes is now 5 and hence 4 legs are required to complete 



a route. Matrix is used to force the algorithm to choose leg one 
with an arc leading from node 1, to one of the other original nodes. 



nodes 2, 3, or 4 , since all other arc choices on the first leg have 



prohibitive costs associated with them. Matrix M^ is a dummy leg 
which is used to form a closed-loop. The only entries of significance 



in M^ are those in the last column, column 5. These m.^ values 
4 i5 

represent the costs of going from any node to node 1, since leg one 



began with an arc leading from node 1. Since the only "acceptable" 



values are m^^ = 5, m^^ = 2, and m^^ = 1 as m^^ = 00 for i = 1 and 5, 
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the optimal route will be forced to close on node 1 as desired. Matrices 

and are identical and are designed to prevent any arc from 

originating at node 1 or node 5 and to prevent any arc from terminating 

at node 1 or node 5, and therefore rows 1 and 5 and columns 1 and 5 

have infinite values. Note that the dotted lines in M and contain 

2 3 

the original matrix less row 1 and column 1, as illustrated by the 
dotted lines in the original matrix. 

The general pattern which emerges is that the matrix for leg 1 
would contain all infinite values except for those arcs leading from 
node 1 to all the other original nodes. The last matrix would contain 
all infinite values except for the last column which would be the same 
as the first column of the original matrix with the infinite value below 
it. The intermediate matrices would be the same as the original matrix 
less row 1 and column 1 with an entire border of infinite values added 
to them . 

With the above modifications, the following traveling salesman 
problems were solved as though they were sequence-dependent routing 
problems. (Only the original matrix is given.) 

Problem No. 4 [Ref. 6] 

1 2 3 4 5 

1 X 5 gTo T 

2 5 X 5 12 12 

3 6 5 X 8 10 

4 10 12 8x6 

5 8 12 10 6 X 

Problem No. 4: Initial Cost Matrix 



Figure 7 
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Problem No. 5 [Ref. 6] 



1 2 3 4 5 6 
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Problem No 


. 5: 


Initial Cost Matrix 


Figure 8 



Problem No. 6 [Ref. 


5] 




1 2 3 4 5 6 
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21 16 25 X 18 18 


5 
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23 5 5 9 5 X 


Problem No 


. 6: Initial Cost Matrix 


Figure 9 
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Problem No . 
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Initial Cost Matrix 



Figure 10 
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Problem No- 8 [Ref. 2] 
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Problem No. 8: Initial Cost Matrix 
Figure 11 



Problem No. 9 [Ref. 3] 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 


X 


24 


18 


22 


31 


19 


33 


25 


30 


26 


2 


15 


X 


19 


27 


26 


32 


25 


31 


28 


18 


3 


22 


23 


X 


23 


16 


29 


27 


18 


16 


27 


4 


24 


31 


18 


X 


19 


13 


28 


9 


19 


27 


5 


23 


18 


34 


20 


X 


31 


24 


15 


25 


8 


6 


24 


12 


17 


15 


10 


X 


11 


16 


21 


31 


7 


28 


15 


27 


35 


19 


18 


X 


21 


21 


19 


8 


13 


24 


18 


13 


13 


22 


25 


X 


29 


24 


9 


17 


21 


18 


24 


27 


24 


34 


31 


X 


18 


10 


18 


19 


29 


16 


23 


17 


18 


31 


23 


X 



Problem No. 9: Initial Cost Matrix 



Figure 12 
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Problem No. 10 [Ref. 7] 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 


X 


28 


57 


72 


81 


85 


80 


113 


89 


80 


2 


28 


X 


28 


45 


54 


57 


63 


85 


63 


63 


3 


57 


28 


X 


20 


30 


28 


57 


57 


40 


57 


4 


72 


45 


20 


X 


10 


20 


72 


45 


20 


45 


5 


81 


54 


30 


10 


X 


22 


81 


41 


10 


41 


6 


85 


57 


28 


20 


22 


X 


63 


28 


28 


63 


7 


80 


63 


57 


72 


81 


63 


X 


80 


89 


113 


8 


113 


85 


57 


45 


41 


28 


80 


X 


40 


80 


9 


89 


63 


40 


20 


10 


28 


89 


40 


X 


40 


10 


80 


63 


57 


45 


41 


63 


113 


80 


40 


X 



Problem No. 10: Initial Cost Matrix 



Figure 13 



Problem No. 1 1 [Ref. 2] 
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Optimal Route "Cost" = 246 

Problem No. 11: Initial Cost Matrix 



Figure 14 
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VI. COMPUTATIONAL RESULTS 



Table I presents summary statistics for the eleven test problems 
considered. Optimal routes obtained verified the known results which 
are in the respective references from which the problems were taken, 
with the exception of test problem Number 7. Reference 2 indicates 
that the optimal route for this problem is as specified in the notes for 
Table I with an optimal route cost of 33. This program obtained the 
optimal route indicated in the table with a route cost of 28 which is 5 
cost units superior to the previous known result. 

The type of problem is either sequence-dependent (SD) or 
traveling salesman problem (TSP) as discussed in Section V. The 
number of complete tours obtained by the program is significant in that 
after the first tour is obtained by branching only to the right, a suc- 
ceeding tour found must have a cost equal to or less than the best tour 
located so far in the computational procedure. It is somewhat repre- 
sentative of the "speed" of convergence towards the optimal solution. 
The number of iterations required is actually the number necessary to 
verify that the best route found by the program is the optimal route. 

The iteration number on which the optimal route is located is in general, 
far lower than the total number of iterations required for verification 
(note test problems Numbers 9 and 10). 

Test problems Numbers 4,5, and 10 have alternate routes 
indicated, but these routes are mirror images of one another and hence 
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are identical. This fact is due to the symmetric nature of the input 
matrices combined with the fact that the matrices are the same for all 
legs of the route excluding the dummy legs. It would be desirable to 
eliminate consideration of any future route which would be an image of 
a route previously located, but this feature was not incorporated into 
the program. For the sequence-dependent symmetric case, test problem 
Number 3, only a single route is found as the matrices for each leg are 
symmetric, but different for each leg. 

Since only a few cases were presented, it would be difficult to 
attempt to draw any conclusions with respect to expected time required 
for. solution of a problem of given size. However, it was observed 
that the time required for a solution rises rapidly as the number of nodes 
increases as discussed in Reference 2. The computer storage require- 
ment for test problems 1 through 9 was a moderate 154K, but problem 
Number 10 required 392K. Test problem Number 11 was run for approxi- 
mately 300 minutes and 25,000 iterations which required 546K bytes 
of storage and the optimal solution was never reached. This matrix is 
symmetric and hence the number of iterations could be reduced by 
taking this fact into account. 

Test problems 7 and 8 terminated in just a few iterations but the 
zeroes in the matrix were placed somewhat strategically. In problems 
9 and 10, the entries in the matrix are nearly all two digit numbers 
which are close to each other in magnitude and hence there is no clear- 
cut minimum route as in problems 7 and 8, and thus the number of 
iterations runs up into the thousands . 
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VII. CONCLUSIONS 



The bra nch-and “bound algorithm and the computer program 
presented can successfully find the optimal route for a variety of 
sequence-dependent routing problems when the matrices of all possible 
costs for each leg of the route are known. 

Although it is admitted that the computer program, as written 
in FORTRAN IV, may not be the most efficient for large-scale problems 
due to storage requirements and processing time, it does provide a 
basis for further programming effort using this algorithm. Although no 
attempt was made to delete nodes from storage once the node bound 
was observed to be above the current least upper bound on a complete 
route, larger scale problems would demand such reduction. 

In the case of symmetric matrices, a programming method must 
be devised to delete consideration of arc (j,i) when arc (i,j) has been 
committed to a leg as this just leads to excessive computation and 
excessive iterations. It is recommended that a lower-level language, 
such as Assembly Language, be utilized to improve efficiency with 
respect to both time and storage requirements since the algorithm 
deals primarily with integer arithmetic operations . 
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APPENDIX A 




COMPLETE COMPUTER SOLUTION TREE FOR TEST PROBLEM NUMBER 2 




u 



APPENDIX B 



CASE NO. 2 

EXAMPLE 2 FRCM REFERENCE 1 
NUMBER OF NODES = 6 

NUMBER OF LEGS = 5 

TYPE problem: SEQUENCE-DEPENDENT 





MATRIX 1 


Ml 






MATRIX 


M2 








MATRIX M3 




40 2A 32 
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12 
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30 18 24 21 9 


36*=!'* 20 26 
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32 
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27 


*** 15 27 3 24 
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32*** 8 
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8** 
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18 


24*** 6 12 12 
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20 20*** 
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16 
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5 
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15 15*** 18 12 
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32 12 8*** 
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2*** 
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24 9 6*** 6 


16 


24 16 2C 
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6 
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5 


3*** 


12 


18 12 15 9*** 





MATRIX 1 


M4 






MATRIX M5 






20 12 16 


14 


6 




50 30 40 35 


15 


18* 


** 10 18 


2 


16 


45=^ 


** 25 45 5 


40 


12 


16*** 4 


8 


8 


30 


40*** 10 20 


20 


6 


10 10*** 


12 


8 


15 


25 25*** 30 


20 


4 


16 6 4*** 


4 


10 


40 15 10*** 


10 


8 


12 8 10 


6* 




20 


30 20 25 15* 





(SEE SOLUTION ON NEXT PAGE) 
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FEASIBLE TCUR NO, 2 IS DECLARED OPTIMAL 
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:{< :!c4c 5}c J}c :{!« ^ 3!< « ajs 5}r :{s :{c :{! ^ :«r :{t :{c 5«c 3^ 

NUMBER OF ALTERNATE OPTIMAL TOURS = 0 

NUMBER OF ITERATIONS REQUIRED = 22 

TIME TO COMPUTE SOLUTION^ 0.738816 SECONDS 



ITERATION INFORMATION 







YTABLE 


DE 


FROM 


WY 


1 


0 


31 


2 


1 


37 


3 


2 


37 


4 


3 


37 


5 


“1 


26 


6 


5 


33 


7 


6 


33 


8 


7 


33 


9 


-5 


30 


10 


9 


37 


11 


“6 


34 


12 


-9 


37 


13 


-12 


41 


14 


-13 


31 


15 


14 


34 


16 


-14 


31 


17 


16 


47 


18 


-15 


35 


19 


-16 


41 


20 


-18 


51 


21 


-19 


32 


22 


21 


59 



YBARTABLE 



10 

2 

3 
6 

4 
2 
3 
6 

5 
2 
3 
3 
2 
2 

3 

4 

5 

3 

4 

5 

4 

5 
3 



JO 

5 

4 
2 
1 

5 

4 

3 
1 

5 

4 

4 

5 

5 
4 
1 
1 
4 

6 
4 
2 
6 
4 



KO TERM 



WY8AR 

20 

35 

45 

32037 

24 

31 
37 

32033 

27 

35 

34 

29 

30 
30 

32 
30 
41 

33 
30 

34 
34 
37 



TERM 
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I 



CASE NO. 6 

EXAMPLE FROM ARTICLE BY LITTLE AND OTHERS IN OR JOURNAL 1963 
NUMBER OF NOCES = 7 

NUMBER OF LEGS = 6 

TYPE PROBLEM: TRAVELING S/^LESMAN CLOSED-LOOP 

MATRICES ARE SAME FOR LEGS 2 THRU (N-2) AND APE AS FOLLOWS: 

I/J= 12 3 4 5 6 7 



1 9999 9999 9999 9999 9999 999^ 9999 

2 9999 9999 16 1 30 25 9999 

3 9999 13 9999 35 5 0 9999 

4 9999 16 25 9999 18 18 9999 

5 9999 46 27 48 9999 5 9999 

6 9999 5595 9999 9999 

7 9999 9999 9999 9999 9999 9999 9999 



FEASIBLE TOUR NO. 3 IS 



LEG FROM TO 

1 1 4 

2 4 3 

3 3 5 

4 5 6 

5 6 2 

6 2 1 



OPTIMAL FOUTE COST = 



DECLARED OPTIMAL 
CCST 
16 
25 
5 
5 
5 
7 
63 



^ ^ 35: ^ sic 3^ ;lc s|c :{t :{t :{s 



NUMBER OF ALTERNATE OPTIMAL TOURS = 0 

NUMBER OF ITERATIONS REQUIRED = 46 

TIME TO COMPUTE SOLUTION= 2.023424 SECONDS 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

^2 

43 

44 

45 

46 



ITERATION INFORMATION 



FROM 


YTABLE 

WY 


10 


JO 


KO 


TERM 


YBAf 

WYBAR 


0 


39 


1 


4 


1 


0 


33 


1 


83 


3 


6 


3 


0 


44 


2 


83 


4 


3 


2 


0 


32083 


3 


88 


2 


5 


5 


0 


99 


4 


88 


6 


2 


4 


0 


32088 


-1 


38 


2 


7 


6 


0 


38 


6 


52 


6 


2 


5 


0 


46 


7 


94 


1 


5 


1 


0 


65 


-2 


60 


3 


6 


4 


0 


49 


9 


65 


2 


3 


3 


0 


71 


10 


65 


4 


2 


2 


0 


32065 


11 


65 


6 


5 


5 


0 


10059 


-6 


43 


5 


7 


6 


0 


46 


13 


50 


1 


6 


1 


0 


44 


14 


65 


3 


5 


5 


0 


63 


15 


74 


2 


4 


3 


0 


80 


-7 


89 


3 


6 


2 


0 


51 


-8 


83 


1 


3 


1 


0 


10021 


-9 


101 


3 


6 


5 


0 


54 


-13 


54 


3 


7 


6 


0 


49 


20 


59 


6 


3 


5 


0 


65 


21 


76 


2 


4 


2 


0 


74 


-14 


45 


1 


2 


1 


0 


60 


23 


63 


2 


4 


2 


0 


60 


24 


68 


4 


6 


3 


0 


70 


-15 


32069 


2 


4 


3 


0 


75 


-17 


83 


3 


6 


3 


0 


56 


-19 


56 


2 


7 


6 


0 


59 


28 


56 


6 


2 


5 


0 


64 


29 


63 


5 


6 


4 


0 


32051 


30 


63 


4 


3 


2 


0 


32063 


-20 


58 


1 


6 


1 


0 


50 


32 


66 


A 


7 


6 


0 


10033 


-23 


62 


1 


3 


1 


0 


10016 


34 


79 


3 


6 


2 


0 


75 


-24 


69 


2 


3 


2 


0 


69 


-27 


32080 


3 


6 


4 


0 


61 


-28 


59 


5 


7 


6 


0 


67 


38 


70 


A 


2 


2 


0 


61 


-32 


51 


1 


2 


1 


0 


53 


40 


69 


2 


4 


2 


0 


66 


-37 


61 


1 


6 


1 


0 


65 


42 


107 


3 


5 


3 


0 


74 


-39 


69 


4 


6 


2 


0 


68 


-40 


58 


1 


5 


1 


0 


66 


45 


66 


5 


6 


2 


0 


80 



TERM 
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APPENDIX C 



COMPOSITION OF COMPUTER CARD DECK 




Variable 

Format 



Variable 

Format 



Variable 

Format 



DATA DECK MAKEUP 



4 

1 


xxxx 

NCASE 

14 


1 


4 6 12 

xxkx!>!^:^ 


13 80 

1 f 


N ALIKE ITS 
14 12 16 


V 

TITLE (I) 
17A4 


4 8 

1 1 


4N 


xxxx xxxx . . . xxxx \ 

M(1,1,1)^M(1,2,1)^. .. 

2014 



Card Type 1 : 
First card in 
Data Deck 



Card Type 2 : 

Second card in Data 
Deck and first card 
of each succeeding 
case 



Card Type 3 : 

Third card on until all matrices 
have been defined. N(N-l) 
cards for each case; input 
matrix for leg 1 first, then 
leg 2 , etc . , row by row. 



Note: All values are integers and must be right- justified in format 
field specified. 
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APPENDIX D 



DESCRIPTION OF VARIABLES USED IN PROGRAM 



PROGRAM CONSISTS OF ALL INTEGER VARIABLES EXCEPT TIMEX 
WHICH IS USED IN CONJUNCTION WITH THE TIMING ROUTINE 

PROGRAM LIMITATIONS : 20 NODES TOTAL INCLUDING THE 

AUGMENTED NODE; NUMBER OF ITERATIONS IS LIMITED BY THE 
NUMBER SPECIFIED BY THE FIRST INDEX OF VARIABLE MATRICES 
YTAB AND YBTAB IN THE DIMENSION STATEMENTS AND MUST BE 
EQUAL TO OR LARGER THAN THE LARGEST VALUE OF THE VARIABLE 
-ITS FOR ANY DATA SET IN THE DATA DECK 



INPUT DECK REQUIREMENTS : CC = CARD COLUMN 

CARD l: F0RMAT(I4) CC 1-4 

NCASE=NUMBER OF CASES TO BE PROCESSED ON THIS RUN 

CARD 2: FORMAT ( 1 4 , I 2 , I 6 , 1 7 A4 ) 

N=NUMBER OF NODES (FORMAT 14) CC 1-4 

ALIKE=1 IF ENTRIES IN MATRIX ARE SAME FOR EACH LEG 
=0 IF ENTRIES IN MATRIX ARE DIFFERENT FOR EACH 
LEG 

(FORMAT 12) CC 5,6 

ITS = MAXIMUM NUMBER OF ITERATIONS DES I RED( FORMAT 16) 
(FORMAT 16) CC 7-12 

TITLE! I) = A HEADING FOR EACH INDIVIDUAL PROBLEM 
(FORMAT 17A4) CC 13-80 



CARD 3 THRU END OF DATA SET: F0RMAT(20I4) 

M(I,J,K) = WORKING SET OF MATRICES: MATRICES ARE 

LOADED ONE ROW AT A TIME 

(MATRIX FOR FIRST LEG, SECOND LEG, ETC,) 

ALL ABOVE DATA IN I FORMAT MUST BE RIGHT JUSTIFIED 
IN FIELD SPECIFIED 

NCASE = NUMBER OF CASES TO BE PROCESSED ON ONE COMPUTER 
RUN 

N = NUMBER OF NODES 

ITS = MAXIMUM NUMBER OF ITERATIONS DESIRED 

ALIKE = 0 IF PROBLEM IS SEQUENCE DEPENDENT TYPE PROBLEM 
ALIKE = 1 IF TRAVELING SALESMAN TYPE PROBLEM 

L = N-1 = NUMBER OF LEGS FOR ROUTE BETWEEN N NODES 

LEGREQ = N-2 = NUMBER OF LEGS REQUIRED TO DETERMINE ROUTE 

COST(K) = COST OF GOING FROM NODE FM(K) TO NODE TO ( K ) 

ON K-TH LEG OF ROUTE 

TCOST(TOUR) = TOTAL COST OF TOUR NUMBER (TOUR) 
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i 



BESKK, J) 



= MATRIX CONTAINING BEST TOUR AT ANY STAGE OF 
SOLUTION AFTER FIRST SUBOPTIMAL TOUR HAS BEEN 
FCUND( K=LEG OF ROUTE) 

J=1 IS LEG OF ROUTE 
=2 IS NODE FROM WHICH LEG K BEGINS 
=3 IS NODE WHICH ENDS LEG K 
=4 IS THE COST TO GO FROM J=2 TO J=3 



BEST(N,1) = NUMBER OF LEAST COST TOUR DETERMINED AT ANY 
POINT AFTER FIRST TOUR IS LOCATED 

BEST(N,2) = CCST OF TOUR NUMBER BEST(N,1) 

LEGCOM(K) = K IF LEG COMMITTED 

= 0 IF LEG UNCOMMITTED 

FM(K) = NODE OF DEPARTURE ON K-TH LEG OF ROUTE 

TO(K) = NODE CF ARRIVAL ON K-TH LEG OF ROUTE 

ARCCOM (I,J) = 100 IF NEITHER I NOR J ARE ON A COMMITTED 
LEG 

=-99 IF EITHER NODE I OR NODE J IS ON A 
COMMITTED LEG 

STEP = STEP NUMBER OF ALGORITHM 
ITER = ITERATION NUMBER 

TOUR = NUMBER CF TOUR FOUND BY ALGORITHM, EACH TOUR 

HAVING A COST WHICH IS LESS THAN OR EQUAL TO THE 
PRECEEDING TOUR 

M(I,J,K) = WORKING SET OF MATRICES 

= COST (OR OTHER VARIABLE TO BE MINIMIZED) OF 
GOING FROM NODE I TO NODE J ON K-TH LEG OF ROUTE 

A(I,J,K) = ORIGINAL M(I,J,K) = PERMANENT FILE OF ALL 
INPUT MATRICES 

MINEL(K) = MINIMUM ELEMENT IN MATRIX K WHEN LEG K IS 
UNCOMMITTED EXCLUDING ROWS AND/OR COLUMNS 
ASSOCIATED WITH NODES ON COMMITTED LEGS 

KIN(K) = CURRENT MINIMUM ELEMENT IN MATRIX K WHEN LEG K 

IS UNCOMMITTED (USED DURING SEARCH FOR MINEL(K)) 

IK(K) = ROW CONTAINING MINEL(K) 

JK(K) = COLUMN CONTAINING MINEL(K) 

THETA = MAXIMUM OF THE SECOND SMALLEST ELEMENTS IN ALL 
RESTRICTED MATRICES FOR UNCOMMITTED LEGS 

MAXEL = CURRENT THETA IN THE DO-LOOP 

MAXLEG = LEG FROM WHICH THETA WAS OBTAINED 

10 = IK(MAXLEG) 

JO = JK(MAXLEG) 

LEG = NCOM = CURRENT NUMBER OF LEGS COMMITTED 

WX=THE LOWER BCUND LABEL ATTACHED TO THE TREE FOR NODE X 

WY=THE LOWER BCUND LABEL ATTACHED TO THE Y NODE OF TREE 

WYBAR = THE LOWER BOUND LABEL ATTACHED TO THE YBAR NODE 
OF THE TREE 
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20 = A LARGE NUMBER ORIGINALLY AND REMAINS AN UPPER BOUND 
ON THE OBJECTIVE FUNCTION 

X(K) = AN ARRAY USED FOR DETERMINING THE FINAL LEG OF 
THE ROUTE AND NODES ON THIS LEG 



INDEX = NODE NUMBER FROM WHICH TO BRANCH 

IS POSITIVE IF BRANCH IS TO BE FROM A Y NODE: 

IS NEGATIVE IF BRANCH IS TO BE FROM A YBAR NODE 

YTABdtJ) = A MATRIX CONTAINING INFORMATION ABOUT Y NODES 
(I) IN COLUMN J WHERE I = ITERATION WHICH 
GENERATED THE NODE 
J = 1 IS THE NODE NUMBER 

= 2 IS THE NODE FROM WHICH BRANCH WAS MADE 
= 3 IS LOWER BOUND LABEL ON NODE Y 
= A IS THE NODE OF DEPARTURE 

= 5 IS THE NODE OF ARRIVAL 

= 6 IS THE LEG OF ROUTE 

= 7 IS ZERO WHEN THE NODE IS NOT TERMINAL 

ONE WHEN THE NODE IS A TERMINAL NODE 

YBTAB(I,J) = A MATRIX CONTAINING INFORMATION ABOUT YBAR 
NODES (I) IN COLUMN J WHERE I = ITERATION 
WHICH GENERATED THE NODE AND IS FOUND IN 
THE MATRIX YTAB{ 1,1) 

J = 1 IS THE LOWER BOUND LABEL ON YBAR NODE 
J = 2 (SAME AS FOR YTABU,7) ) 

(NOTE: THE NODE NUMBER IS THE NEGATIVE OF THE 

CORRESPONDING Y NODE FOR THE SAME ITERATION, I.) 



BB, G, FROM, NCOM : ALL ARE VARIABLES USED IN RECON- 
STRUCTING MATRICES WHEN NODE FROM WHICH BRANCH IS TO 
OCCUR IS NET NODE OF PREVIOUS STEP 6 

KEY : IS A VARIABLE USED TO CONTROL FLOW OF PROGRAMMING 

THROUGH ALGORITHM TO AVOID ADDITIONAL DUPLICATE 
CODE WHICH WOULD BE REQUIRED 

GOLF, DEL : VARIABLES USED TO DENOTE CURRENT ROW OR 
COLUMN OF ARCCOM MATRIX WHICH MAY BE USED FOR 
NOT CONSIDERING CERTAIN ELEMENTS OF THE MATRICES 

FUNCTION SUBPROGRAMS USED : 

MINO - FINDS MINIMUM OF 2 OR MORE INTEGER^A ARGUMENTS 
AND ASSIGNS A FUNCTIONAL INTEGER VALUE 

MAXO - FINDS MAXIMUM OF 2 OR MORE INTEGER*4 ARGUMENTS 
AND ASSIGNS A FUNCTIONAL INTEGER VALUE 
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APPENDIX E 

COMPUTER PROGRAM LISTING 



:(e 3{t 5[< 5jc :i!t :{! :<t ^ :?« :5c s^e ^ -c J}c ^ 3|< :{e :4s j}£ 



4s 4s 

* A COMPUTER PROGRAM FOR THE SOLUTION 

OF SEQUENCE-DEPENDENT ROUTING PROBLEMS =4= 
’f' USING A BRANCH-AND-BOUND ALGORITHM * 

4s 4c 

« :{c 



4s4s4s4«^s;<:<<5<c3^4s 4s 4s 4s 4s 4c 4s 4s 4: 4s 4c 4s 4s 4: 4c 4c 4s 4c 4s 4s 4s 4s 4s 4s 4c 4s 4s 4s 4s 4s 4c 4s 4s 4s 4: 4s 4c 



Q 4t5?ci!c 4=5^ Jk^«5i'4:5;cs5:2{c4c5!:4«5!c=<c5!«^^5!c*’«c5!<^5{c5!<5{c3;«4<5f<5r5!<^>;c***>{c5{c5l‘i!c=(«os4s4s4s4s4c4c4s4s4s4s>Jr4c4c 



IMPLICIT INTEGER^2 ( A-Z) 

REAL*4 TIMEX 

INTEGER’^A M( 20, 20, 19) ,THETA, MAXEL , MINEL ( 20) ,MIN( 20) , 
1TITLE(17),ZC,WY 

DIMENSION A (20,20, 19) , COST (20) ,LEGC0M(20 ) ,TC0ST(30) 
1FM(20),T0(2C),X(20) , ARCCOM ( 20 , 20 ) , BEST ( 20,4) 
2YTAB(2500,7) , YBT AB ( 2500 , 2 ) , I K ( 20 ) , JK ( 20 ) 

COMMON TIMEX,M,THETA,MAXEL,MINEL,MIN,TITLE,ZO,WY,A, 
ICOST, LEGCCM,FM,TO, X , ARCCOM , B E ST , YT AB , Y BTAB , IK, JK, TOUR, 
2AA,N,L, ALIKE, LEGREQ, STEP, I TER, DEL, WX, MAX LEG, LEG, WYBAR, 
3I0,J0,KC,ITS, TCOST, INDEX 

1 FORMAT (14) 

2 FORMAT ( 14, 12,16, 17A4) 



READ(5,1) NCASE 
DO 2000 AA = 1, NCASE 
C READ INPUT PARAMETERS 

READ (5, 2 ) N, ALIKE, ITS, (TITLE ( I ), 1=1, 17 ) 



CALL INPUT 

CALL ITERTE (C2000) 

C OPTIMALITY REACHED : PROCESS A NEW CASE OR TERMINATE 

CALL SOLN 

2000 CONTINUE 
STOP 
END 
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SUBROUTINE INPUT 



IMPLICIT INTEGER^2(A-Z) 

REAL=!^4 TIMEX 

INTEGER’S'A M( 20,20, 19) , THETA , MAXEL , MI NEL ( 20 ) ,MIN(20) , 
1TITLE( 17),ZC,WY 

DIMENSION A (20,20, 19) , COST (20) ,LEGC0M(20) ,TC0ST(30) 
1FM(20) ,T0( 20 ,X(20) ,ARCC0M(20,20) ,BEST(20,4) 
2YTAB(2500, 7 ) , YBTAB ( 2500 , 2 ) ,IK(20),JK(20) 

CCMMCN TIMEX, M, THETA, MAX EL, MINEL, MIN, TITLE, ZO,WY, A, 
ICOST ,LEGCOM, EM, TO, X,ARCCOM, BEST, YTAB, YBTAB, IK, JK, TOUR, 
2AA,N,L, ALIKE, LEGREQ, STEP, I T ER , DEL , WX , MAX L EG , LEG , W YB AR , 
3I0,J0,KC,ITS, TCOST, INDEX 

DATA HEAD/1 ,2 ,3,4, 5,6, 7,8,9, 10, 11, 12 , 13 , 14 , 1 5 , 16, 1 7 , 
118, 19,20,21,22,23,24/ 



2 FORMAT (2014) 

4 FORMAT (IHl) 

5 FORMAT( • 0* , lOX, 'CASE NO. * , 13 ,// , 11 X, 17A4, //I IX , * NUMBE* 
1,*R OF NODES = * ,12,//, IIX, 'NUMBER OF LEGS = *,I2) 

6 FORMAT (• O' , lOX, * TYPE PROBLEM: SEQUENCE-DEPENDENT') 

7 FORMAT (' C , lOX, ' TYPE PROBLEM: TRAVELING SALESMAN', 

1' CLOSLD-LCCP' ) 

8 FORMAT( 'O' , 15X, 'MATRIX M 1 ' , 1 4X , ' M ATR I X M2',14X, 
I'MATRIX M3 ' ,14X, ' MATRI X M4 ' , 14X , * M ATRI X M5') 

9 FORMAT (/ ,11X,6I3,5X,6I3,5X,6I3,5X,6I3,5X,6I3) 

10 FORMATCO' ,20X, 'MATRIX M 1 ' , 20 X , ' M ATR 1 X M 2',20X, 
I'MATRIX M 3 ', 20X, • MATRIX M 4') 

11 F0RiMAT(/,10X,5I5,5X,5I5,5X,5I5,5X, 515) 

33 FORMAT( IIX, 12,2/ ,2215) 

37 FORMAT ( 'O' , lOX, 'MATRICES ARE SAME FOR LEGS 2 THRU ', 
l'(N-2) AND ARE AS FOLL OWS : ' / / , 1 1 X , ' I / J = ' , 2 2 I 5 ) 

39 FORMAT! '0* ) 

C INITIALIZATION 

DO 101 1=1, N 
COST( I ) = 0 
LEGCCM( I ) = 0 
FM( I ) = 0 
T0( I ) = 0 
DO 101 J=1,N 
101 ARCCGM( I , J ) = 100 
DEL = 1 
ITER = 0 
TOUR = 0 
LEG = 0 
L = N-1 

LEGREQ = N - 2 

C NOTE: VALUES OF 32000 IN THE PROGRAM ARE USED TO 
C INDICATE INFINITE VALUES 

ZO = 32000 

STEP = 1 

C READ INPUT MATRICES FOR ALL L LEGS 

IF (ALIKE. EQ.l) GO TO 104 
DO 103 K=1,L 
DO 103 1=1, N 

103 READ(5,2) ( M ( I , J , K ) , J= 1 , N ) 

GO TO 110 
1C4 DO 1C5 K=l,2 
DO 105 1=1, N 

105 READ(5,2) ( M ( I , J , K ) , J= 1 , N ) 
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DO 106 K=3,LEGREQ 
DO 106 1 = 1, f\ . 

DO 106 J=1,N 

106 M( I , J,K) = N( I , J,K-1) 

DO 1C7 1 = 1, N 

107 READ(5,2) ( N ( I , J , L ) , J= 1 , N ) 

C WRITE OUT INPUT 

110 WRITE(6,4) 

IF( ALIKE. EQ.O) WRITE(6,6) 

IF( ALIKE. EQ.l ) WRITE(6,7) 

WRITE! 6, 5) AA, (TITLE ( I ) , 1 = 1, 17 ) ,N, L 
IF(N.NE.5.AND.N.NE.6.0R.ALIKE.EQ.l) GO TO 113 
IF(N.E0.5) WRITE(6,10) 

IF(N.EQ.6) WRITE(6,8) 

DO 111 I=1,N 

IF(N.EO. 5) WR ITE(6, 11 ) ( ( M ( I , J , K ) , J= 1 , N ) , K = 1 , L ) 

111 IF(N.EQ.6) WRITE (6,9) ( ( M ( I , J , K ) , J = 1 , N ) , K= 1 , L ) 

GO TO 131 

113 WRITE(6,37) (HEAD! I ) , I = 1,N) 

WRITE(6,39) 

DO 117 1=1, N 

117 WRITE(6,33) I , ( ( M ( I , J , K) , J = 1 ,N ) , K=2 , 2 ) 

131 CONTINUE 

C START TIMER 

CALL TIMEIT (0, TIMEX) 

C CREATE COPY OF ORIGINAL DATA IN MATRIX A 

112 DO 132 K=1,L 
DO 132 1=1, N 
DO 132 J=1,N 

132 A(I, J,K) = P( I,J,K) 

C STEP TWO 

C FIND MINIMUM ELEMENT IN EACH MATRIX 

DO 200 K=1,L 

IK(K) = 0 

JK(K) = 0 

MIN(K) = 32C00 

DO 190 1=1, N 

DO 190 J =1,N 

IF( I .EQ. J) GO TO l<^0 

MIlMEL(K) = MINO ( M I N ( K ) , M ( I , J , K ) ) 

116 IF (MINEL(K).GE.MIN(K) ) GO TO 190 
IK(K) = I 
JK(K) = J 
MIN(K) = MINEL(K) 

190 CONTINUE 
2C0 CONTINUE 

C REDUCE MATRICES 

DO 212 K=1,L 

DO 212 1=1, N 

DO 212 J=1,N 

IF(I.EQ.J) GO TO 212 

M( I , J,K) = M( I , J,K) - MINEL(K) 

212 CONTINUE 

C LABEL NODES 

WX = 0 

DO 230 K=1,L 
230 WX = WX + MINEL(K) 

RETURN 

END 
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SUBROUTINE ITERTE («) 



IMPLICIT INTEGER=«'2(A-Z) 

REAL’S^A TIMEX 

INTEGER^A M(20t20, 19) » THETA t MAXEL , MI NEL ( 20 ) ,MIN(20) , 
1TITLE( 17) ,ZC,WY 

DIMENSION A ( 20 t 20, 19) t CO ST (20) ,LEGC0M(20) ,TC0ST(30) 
1FM( 20) ,T0( 20 ,X( 20) , ARCCOM( 20 ,20) ,8EST( 2 0,4) 

2YT AR ( 2 500,7) ,YBTAB( 2500,2) , I K ( 20 ) , JK ( 20 ) 

COMMON T I M F X , M , T HE T A , M AX EL , M I NE L , M I N , T I T L E , Z 0 , WY , A , 
ICOST, LEG COM, F M, T 0, X,ARCC0M, BEST, YTAB,YBTAB, IK, JK, TOUR, 
2AA,N,L , ALI KE, LEGREQ,STEP, I TER , DEL , WX , MAX L EG , LEG , WYB AR , 
310, JO, KO, ITS, TCO ST, INDEX 

12 FORMAT! * 0* , lOX, * MAXIMUM ARRAY STORAGE SPACE EXCEEDED:* 
1,* ITER=*,T6,* AND MATRICES YTAB AND YBTAB ARE DIMEN*, 

2*SI0NED FOR*, 17* I T FR A T I ON S . * , / / , 1 IX , * NUMB ER OF *, 
3*T0URS OBTAINED =*,I3,». BEST VALUE SO FAR IS *,I5, 

4* FOR TOUR NUMBER* , 13, *.* ) 

15 FORMAT! *0* , lOX, * THE INFORMATION OBTAIN BY THE PROGRAM* 
1,* SO FAR IS PRINTED OUT BELOW, SOLELY FOR REFERENCE*, 
2*. BEST ROUTE HAS NOT BEEN FOUND.*) 

STEP THREE - ITERATION PROCEDURE 

CONSISTS OF STEPS THREE THRU EIGHT 

LABEL = 1 

5C0 I TER = ITER + 1 

C MAXIMUM ARRAY STORAGE SPACE EXCEEDED: ERROR 

IF! ITER. GT. ITS) CALL TIMEIT (-1, TIMEX) 

IF! ITER. GT. ITS) WR I TE ( 6 , 12 ) I TER , I T S , TOUR , B EST ( N , 2 ) , 
1BEST!N,1 ) 

IFdTER.GT. ITS) GO TO 840 

MAXEL = -10 
MAXLEG = 0 
DO 512 K=1,L 

IF(K.EQ.LEGCOM!K) ) GO TO 512 

MIN! K) = 32C00 

MINEL!K) = -1 

DO 509 1=1, N 

DO 508 J=1,N 

IF(I.EQ.J) GO TO 508 

IF (K.EQ.L) GO TO 505 

IF!K.NE.LEGC0M!K+1)-1) GO TO 505 

IF!K.EQ.LEGC0M(K-1 )+l) GO TO 510 

IF! J.NE. FM! K+1 ) ) GO TO 508 

IF! ARCCOM! I ,DEL) .LT.lOO) GO TO 508 

GO TO 504 

5C5 IF(K.EQ.l) GO TO 506 

IF!K.NE.LEGC0M!K-1)+1) GO TO 506 
IF! I .NE.TO! K-1 ) ) GO TO 508 
IF!ARCCOM! DEL, J) .LT.lOO) GO TO 508 
GO TO 504 

506 IF!ARCCOM( I , J ).NE. 100) GO TO 508 
IF!ARCCQM! J, I ) .NE. 100) GO TO 508 
504 IF ! I.EQ.IK!K) .AND. J.EQ. JK!K) ) GO TO 508 
MINEL!K) = MINO!MIN!K) ,M!I ,J,K) ) 

IF!MINEL!K) .GE.MIN(K) ) GO TO 508 
MIN! K) = MINEL!K) 

5C8 CONTINUE 
5C9 CONTINUE 

510 IF! MINEL!K) .EQ.-l ) MINEL!K) = 32000 

THETA = MAXC!MAXEL,MINEL(K) ) 

IF!THETA.LE. MAXEL) GO TO 512 
MAXLEG = K 
MAXEL = THETA 
512 CONTINUE 
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C STEP FOUR ITERATION PROCEDURE 

C LABEL YBAR BY \nYBAR 

WYBAR = WX + THETA 
LEGCCM(MAXLEG) = MAXLEG 
FM(MAXLEG) = IK(MAXLEG) 
TO(MAXLEG) = JK(MAXLEG) 

JO = TO(MAXLEG) 

10 = FM(MAXLEG) 

CGST(MAXLEG) = A ( I 0 , JO , M AXL EG ) 



YBTAB( ITER,1) = WYBAR 
YBTAB( ITER, 2) = 1 

YTABdTER,!) = ITER 

IF( ITER.EQ. 1 ) YTABdTER, 2) = 0 

IF( ITER.GT . 1 ) YTAB(ITER,2) = INDEX 

YTABdTER, A) = 10 

YTABdTER, 5) = JO 

YTABdTER, 6) = MAXLEG 

YTABdTER,?) = 1 

IF LEG JUST DETERMINED ALLOWS A ROUTE TO BE SPECIFIED, 
GO TO STEP 7 

IF(LEG.EQ.N-3) YTABdTER, 3) = WX 
IF(LEG.EQ.N-3) GO TO 700 



C STEP FIVE ITERATION PROCEDURE 

C DELETION OF ARCS NOT POSSIBLE 

DO 513 1=1 ,N 
DO 513 J=1,N 
IFd .EQ. J) GO TO 513 

I F( ( I. EQ.IO.OR. I .EQ. JO.OR. J.EQ. lO.OR, J.EQ. JO) .AND. 
1ARCCCM( I , J ) .GT.99) ARCCOM(I,J) = -99 
513 CONTINUE 



GOLF = 1 

515 DEL = GOLF 
DO 516 F=1,L 

IF (FM(F).FO.DEL) GOLF = GOLF + 1 
IF (FM(F).EC.DEL) GO TO 515 
IF (TO(F).EG.DEL) GOLF = GOLF + 1 

516 IF (TO(F).EC.DEL) GO TO 515 

C STEP SIX - - ITERATION PROCEDURE 

C FIND MINIMUM ELEMENT IN EACH MATRIX K WHERE LEG K 
C IS UNCOMMITTED 

DO 630 K=1,L 
KEY = 0 
IK(K) = 0 
JK(K) = 0 

IF (K.EQ.LEGCOM(K) ) GO TO 630 
MINEL(K) = -1 
602 MIN(K) = 32000 
DO 629 1=1, N 
DO 629 J=I,N 
IF( I .EO.J) GO TO 621 

IF ( K.EQ.l . /5ND.LEGCGM( 2) .EQ.O) GO TO 605 
IF ( K.EQ.l ) GO TO 606 

IF (K.EQ.L. /^ND.LEGCOM(L-l) .EQ.O) GO TO 605 
IF (K.EQ.L) GO TO 610 

IF( K.EQ. LEGC0M(K+1 )-l. AND. K. NE . LEGCOM ( K- 1 ) +1 )GOTO 606 
IF( K. EQ. LEGCOM (K-1 ) +1. AND. K.NE .LEGC0M(K+1 ) -1 ) GOTO 610 
IF(K.EQ. LEGCOM (K+1 )-l. AND. K. EQ. LEGCOM (K-1 ) +1 )GOTO 609 
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605 IF ( APCCCM( I, J).LT*100) GO TO 621 
IF (KEY.EQ.C) GO TO 620 

IF ( KEY. EQ. 1) GO TO 622 

606 IF ( J.N£.FM(K+n ) GO TO 621 

IF ( ARCCCMI It DEL ) .LT. 100) GO TO 621 
IF (KEY.EQ.O) GO TO 620 
IF (KEY.EQ.l) GO TO 622 

609 IF( I .EO.TO (K-1 ) . AND.J. EQ.Fn(K+l ) . AND . KEY . EQ. 0 ) G0TG620 
I F( I .EO.TO( K-1 ) . AND. J. EQ.FMI K+1 ) .ANO.KEY.EQ. 1 )G0T0622 
GO TO 621 

610 IF ( I.NF.TO(K-l) ) GO TO 621 

IF ( APCCOM.IOFL, J ) .LT.lOO) GO TO 621 
IF (KFY.EQ. 1) GO TO 622 

620 MINFL(K) = I NO ( N I N ( K ) t M ( I , J t K ) ) 

IF (mNEL(K).GE.MIN(K) ) GO TO 621 
IK(K) = I 

JK(K) = J 
HIN(K) = MINEL(K) 

621 IF (KEY.EQ.l) GO TO 629 
IF ( I.NE.N) GO TO 629 
IF (J.NE.N) GO "^0 629 

C REDUCE MATRICES 

KEY = 1 
GO TO 602 

622 IF(I .EQ. J) GO TO 629 
M(lTJtK) = M(I,J,K) MINEL(K) 

629 CONTINUE 

630 CCNTINUE 

C LABEL Y BY WY 

WY = WX 
DO 635 K=1tL 

IF(K.EO.LEGCGM(K) ) GO TO 635 
WY = WY + MINEL(K) 

635 CONTINUE 

YTAB(ITER,3) = WY 

STEP = 7 - - INCREMENT NUMBER OF LEGS COMMITTED AND 
DETERMINE V^HAT STEP IS NEXT 

700 LEG = LEG + 1 

IF(LEG.NE.N/2) GO TO 720 
CALL CHECK(S720) 

CALL ROUTE (6800,0840) 

720 IF(LEG.GE.LEGREQ.AND.WY.LE.ZO) CALL ROUTE ( 6800 , 6 840 ) 
IF(LEG.GE.LEGREO.AND.WY.GT.ZO) GO TO 799 
IF(LEG.LT.LEGREQ.AND.WY.LE.ZO) GO TO 850 



C STEP 8 - SELECT NODE X FROM WHICH TO BRANCH 

C MODIFICATION TO ORIGINAL ALGORITHM FOR DETERMINING BRANCH 
C POINT - - BRANCH TO THE RIGHT WHENEVER A TOUR IS NOT 
C COMPLETED OR BRANCH FROM THE LOWEST NUMBERED YBAR NODE 
C Vi/HICH IS A TERMINAL NODE, IN THE ORDER GIVEN. 



799 YTAB(ITER,7) = 0 
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C BRANCH TO RIGHT IS EXHAUSTED, THEREFOR SEARCH YBAR NODES 
C FOR A FEASIBLE LABEL (I.E. LABEL. LE.ZO) 

800 DO 830 I=LABEL,ITER 

IF(YBTAB(I,^).EO.O) GO TO B30 
IF(YBTAB(I ,3) .GT.ZO) YBTAB(I,4) = 0 
IF{ YBTABl I , 3) .GT.ZOI GO TO 830 
LABEL = I 

INDEX = YBTABI 1,1) 

CALL SETUP (8500) 

830 CONTINUE 

C OPTIMAL ROUTE HAS BEEN FOUND : TERMINATE 

CALL TIMEIT (-1, TIMEX) 

RETURN 



C . ERROR MESSAGES HAVE BEEN PRODUCED: TERMINATE CASE 

840 WRJTE(6,15) 

RETURN 

C BRANCH TO RIGHT IS NOT EXHAUSTED: THEREFOR BRANCH FROM 
C Y NODE AND MAKE NODE.Y NON-TERMINAL 



850 YTAB( ITER,7) = 0 
INDEX = YTAB( ITER, 
WX = WY 
GO TO 500 
END 



1 ) 
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SUBRCUTINE SETUP (=«'') 



IMPLICIT INTEGER*2 ( A-Z ) 

REAL={=4 TIMEX 

INTEGER^4 M(20t20,19), THET A , MAXEL , MI NEL ( 2 0 ) ,MIN(20) , 
ITITLEI 17) t ZC, WY 

DIME NS ION A (20,2 0, 19) ,COST ( 20) ,LEGCQM(20) , TCOST( 30) 
1FM(20),T0(2C),X(20) , ARCCCM ( 20 , 20 ) ,BEST(20,4) 
2YTAB(25C0,7) ,YBTAB( 2500,2) , I K ( 20 ) , JK { 20 ) 

COMMON TIMEX, M,THET A, MAXEL , MI NEL, MIN, TITLE, ZO,WY, A, 
1C0ST,LEGC0M,FM,T0, X,ARCCOM,BEST, YTAB,Y6TAB, IK, JK, TOUR 
2AA,N,L, ALIKE, LEGREQ, STEP, I TER, DEL, WX, MAX LEG, LEG, WYBAR 
310, JO, KO, ITS, TCO ST, INDEX 

C STEP NINE 

LEG = 0 
G = 0 
NCCM = 0 

IF( INDEX. LT.O ) Y BT AB (- INDEX , 2 ) = 0 

IF( INDEX. GT.O ) Y TAB ( I N DEX , 7 ) = 0 

DO 901 1=1, N 

LEGCOM( I ) = 0 

FM( I ) = 0 

TOC I ) = 0 

901 COSTd ) = 0 
DO 902 1=1 ,N 
DO 902 J=1 ,N 

902 ARCCGM( I , J ) = 100 

C STEP 9 SUBSTEP 1 

C COMPUTE G=SUM A(I,J,K) FOR COMMITTED ARCS AND LEGS 
BB = INDEX 

IF( INDEX. LT.O) FROM = YTAB (-INDEX, 2 ) 

IF( INDEX. GT.O ) FROM = YT AB ( I NDEX , 2 ) 

DO 909 1=1 , ITER 
IF (PB.GT.O) GO TO 903 

IF ( FROM.EC.-l .OR.FROM.EQ.O) GO TO 910 
IF (FROM. LT.O) GO TO 908 
GO TO 904 

C LSED TO START EACK TREE FROM BRANCH NODE AND 
C CONSIDER THAT NODE 

903 FROM = BB 



904 G = G + A(YTAB(FR0M,4) ,YTAB(FR0M,5),YTAB(FR0M,6) ) 
LEGCCMC YTAB(FR0M,6) ) = YTAB(FR0M,6) 

KO = YTABC FROM, 6 ) 

FM(YTAB( FR0M,6) ) = YTAB(FR0M,4) 

TO = YTABC FROM, 4) 

TC(YTAB(FR0M,6) ) = YTAB(FR0M,5) 

JO = YTAB( FROM, 5 ) 

COST(YTAB( FR0M,6) ) = A(IO,JO,KO) 

DO 905 AI=1,N 
DO 905 BJ=1,N 

IF( ( AI . EQ. IC .OR. AI . EQ. JO.OR.BJ .EQ. lO.OR. BJ.EQ. JO) .AND 
1ARCC0M(AI,BJ) .GT.99) ARCCOM ( A I , B J ) =-99 

905 CONTINUE 

NCGM = NCOM + 1 
FROM = YTAB(FR0M,2 ) 

BB =-1000 
GO TO 909 

908 FROM = YTAB (-FROM, 2 ) 

909 CONTINUE 



52 



C STEP 9 SUBSTEP 2 



SETTING UP M(K) 



C 

C 



C 



C 



+ 1 
+ 1 



910 LEG = NCCM 
DO 911 K=1,L 
DO 911 1=1, N 
DO 911 J=1,N 

911 M( I, J,K) = A( I ,J,K) 

STEP 9 SUBSTEF 3 
DELETE ARCS AND LEGS COMMITTED 

GOLF = 1 

912 DEL = GOLF 
DO 913 F=1,L 

IF {FM(F).EG.DEL) GOLF = GOLF 
IF (FM(F).EC.DEL) GO TO ^12 
IF (TQ(F). EC.DEL) GOLF = GOLF 
IF (TO(F).EC.DEL) GO TO 912 

913 CONTINUE 

STEP 9 SUBSTEP A BLOCK PATHS NOT ALLOWED 

IF( INDEX.FQ.-l )M(YTAB( 1 ,4) ,YTAB( 1, 5 ) , YTAR( 1 , ^ 
IF( INDEX. EQ.-l ) GO TO 919 

IF (INDEX. LT.O) M ( YT AB ( - INDEX , 4 ) , YTAB ( - I ND EX , 
1-INDEX, 6) ) = 32000 
IF(INDEX.LT.O) FROM = YT AB ( - I NDE X , 2 ) 

IF( INDEX. GT.O ) FROM = YTAB ( I NDEX , 2 ) 

DO 918 1=1 , ITER 
IF ( FROM. EG. 0) GO TO 919 
IF (FROM. EQ.-l) GO TO 916 
IF (FROM. LT.O) GO TO 917 
FROM = YTAB(FR0M,2 ) 

GO TO 918 

917 M( YTAB (-FROM, 4 ) , YT AB ( -FROM , 5 ) , YT AB ( -FROM , 6 ) ) 
FROM = YTAB (-FROM, 2 ) 

918 CONTINUE 

916 M(YTAB( 1,4) ,YTAB(1,5) ,YTAB( 1,6) ) = 32000 



STEP 9 SUBSTEP 5 



FIND MINIMUM ELEMENT IN EACH 



919 CALL MINELM 



WX = G 

DO 940 K=1,L 

IF ( K.EQ.LEGCOM(K) ) GO TO 940 
WX =WX + MINEL(K) 

940 CONTINUE 
RETURN 1 
END 



) ) = 32000 
5),YTAB( 

= 32000 
MATRIX 
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SUBROUTINE ROUTE (*»-) 



IMPLICIT INTEGER*2(A-Z) 

TIMEX 

INTEGER M(20,20» , THETA , MAXEL » MI NEL ( 20 ) »MIN(20) , 

1TITLE( 17) , ZC,V>IY 

DIMENSION A (20,20, 19) , COST ( 20) ,LEGC0M(20) ,TC0ST(30) 
1FM(20) ,T0( 20) , X(2Q) , ARCC GM ( 2 0 , 20 ) , BEST ( 20,4) 
2YTA8(2SOO,7 ) , YBT AB { 2500 , 2 ) , I K ( 20 ) , JK { 20 ) 

COMMON TIMEX,M,THETA, MAXEL,MINEL,MIN,TITLE, ZO, WY, A, 
1C0ST,LEGC0M,FM,TG, X,ARCCOM, BEST, YTAB, YBT AB, IK, JK, TOUR, 
2AA,N,L,ALIKE, LEGREg,STEP, ITER ,OEL,WX, MAXLEG,LEG,'JYBAR, 
3IO,JO,KO,ITS, TCOST, INDEX 

C STEP 10 

. 20 FORMAT! '0* ,10X,* THE TOUR NUMBER IS EQUAL TO 31, AND*, 
1* THE VARIABLE TCOST IS ONLY DIMENSIONED FOR 30 TOURS* 
2, * : CASE TERMINATED. * ) 

21 F0RMAT(*0* ,5X,*STEP NO.*, 13,* ITER NO.*, 16,*: UPPER*, 
1* BOUND ON VALUE OF OPTIMAL TOUR IS*, 15) 

23 FORMAT ( * 0 * , 5X, * FE ASI BLE TOUR NO.*, 13,* IS AS FOLLOWS: 



ZO = WY 

WRITE! 6,21 ) STEP, ITER, ZO 

C NODE IS MADE NCN-TERMINAL SINCE A ROUTE HAS BEEN 
C COMPLETED AND NO BRANCHING CAN TAKE PLACE. 

YTAB!TTER,7) = 0 
TOUR = TOUR + 1 

C ERROR MESSAGE: THE NUMBER OF TOURS IS. GT. TCOST DIMENSION 

IF(T0UR.GT.30) WRITE!6,20) 

IF(TCUR.GT.30) RETURN 2 

DO 985 1=1, K 
985 X! I ) = 0 

C DETERMINE BY PROCESS OF ELIMINATION AND ORDERING 
C WHAT LEG OF ROUTE IS MISSING AND THUS FORM COMPLETED 
C ROUTE. 

1000 DO 1020 K=1,L 

IF !K.EQ.LEGCOM!K) ) GO TO 1020 
LEGCOM(K) = K 
IF ! K.NE.l ) GO TO 1010 
T0!1) = FM(2) 

DO 1002 1=1, L 

IF (FMUl.NE.O) X!FM!D) = FM! I ) 

1002 CONTINUE 
X!TO(L) ) = TO!L) 

DO 1003 1=1, N 

IF ! X! I ) .NE.O) GO TO 1003 
FM! 1 ) = I 

C0ST!1) = A!FM!1),T0!1),1) 

GO TO 1021 

1003 CONTINUE 
GO TO 1020 

1010 IF (K.EQ.L) GO TO 1011 
FM!K) = TO(K-l) 

TO(K) = FM!K+1) 

COST(K) = A!FM!K) ,TO!K) ,K) 

GO TO 1021 

1011 FM(L) = T0!L-1) 

DO 1012 1=1, L 

IF !TO!I).NE.O) X!TO(D) = TO(I) 

1012 CONTINUE 
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X(FM(1) ) = FM(1) 

DO 1013 1=1, N 

IF ( X( I ) .NE.O) GO TO 1013 

TQ(L) = I 

COST(L) = A{FM(L),TO(L) ,L) 
GO TO 1021 
1013 CONTINUE 
1020 CONTINUE 



1021 TCOST(TOUR) = 0 
DO 1030 K=l,L 

1030 TCOST(TOUR) = TCOST(TOUR) + COST(K) 

C COMPLETE TOUR IS NOW KNOWN. ENTER IT IN MATRIX BEST. 

IF(TOUR.EQ. n GO TO 1040 
IF(TC0ST(TQUR).GE.BEST(N,2M GO TO 1060 
lOAO DO 1050 K=1,L 
BEST (K,l) = K 
BEST (K,2) = FM(K) 

BEST (K,3) = TO(K) 

1050 BEST (K,4) = COST(K) 

BEST(N,1) = TOUR 
BEST(N,2) = TCOST(TOUR) 

1060 RETURN 1 
END 
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SUBROUTINE SGLN 



IMPLICIT INTEGER^2 ( A-Z) 

REAL=f=4 TIMEX 

INTEGER'^4 M (20, 20 1 19) , THETA , MAXEL , MI NEL ( 20 ) ,MIN(20) , 
1TITLE( 17) , ZC, WY 

DIMENSION A (20, 20, 19) , COST (20) ,LEGC0M(20) ,TC0ST(30) 
1FM( 20) , TG( 20 , X( 20) , ARCCOM( 20,20) ,BEST( 20, A) 
2YTAB(2500, 7 ) , YBTAB ( 2500,2) ,IK(20),JK(20) 

COMMON TIMEX,M,THETA,MAXEL, MINEL,MIN,TITLE,Z0,WY, A, 
1C0ST,LEGC0M,HM,T0, X,ARCCOM, BEST, YTAB, YBTAB, IK, JK, TOUR, 
2AA,N,L ,ALIKE,LEGREQ,STEP, ITER,DEL, WX,MAXLEG,LEG,WYBAR, 
3IO,JO,KO,ITS, TCOST, INDEX 

4 FORMAT! IHl) 

5 F0RM/iT(*0* ,///,25X, 'ITERATION INFORMATION') 

15 FORMAT (' O' , lOX 'NUMBER OF ALTERNATE OPTIMAL TOURS =', 

1 14) 

16 FORMAT! 'O' , lOX, ' BEST TOUR SO FAR IS AS FOLLOWS:') 

17 FORMAT! 'O' ,26X, 'ROUTE COST = ',15) 

18 FORMAT ( 'O' , lOX, 'FEASIBLE TOUR NO. ',13,' IS DECLARED', 
1' OPTIMAL') 

19 FORMAT ( • 0 ' , ICX , ' ' ) 

24 FORMAT!' ' , 7X , 31 8, 415, 4X,2 18 ) 

27 FORMAT! ' O' , lOX ,' LEG' ,5X, 'FROM' ,5X, ' TO ' , 5X , » COST * ) 

2 8 FORMAT! 'O' , IIX, 12, 7X, I 2, 6X, 12, 5X, 14) 

31 FORMAT !• O' ,27X, ' YTABLE ' , 27X , ' Y6 ART ABLE' , / / , 12X, 

I'NQDE FRCM WY 10 JO KO TERM',7X, 

2'WYBAR TERM') 

48 FORMAT! ' O' , lOX, ' NUMBER OF ITERATIONS REQUIRED =',I7) 

50 FORMAT! 'O' ,13X, 'OPTIMAL ROUTE COST = ',15) 

55 FORMAT! ' O' , lOX, ' TIME TO COMPUTE SOLUT I 0N= ' ,-6PF 1 5 . 6 , 

1' SECONDS') 

56 FORMAT! ' O' , lOX, ' TI ME USED UP SO FAR = ',-6PF15.6, 

1' SECONDS') 



IF(ALIKE.EQ.l) BEST!L,3) = 1 
IF! ITER. GT • ITS) GO TO 1520 
WRITE! 6, 19) 

WRITE!6,18) BEST!N,1) 

1520 IF! ITER. GT. ITS) WR1TE!6,16) 

WRITE! 6, 27 ) 

DO 1550 K=1,L 

1550 WRITE(6,28) ! B EST ! K , J ) , J= 1 , 4 ) 

IF! ITER. LE. ITS) WRITE (6,50) BE ST (N, 2) 
IF! I TER. GT. ITS) WRITE! 6, 17) BEST!N,2) 
WRITE!6, 19) 

IF! ITER. GT. ITS) GO TO 1575 
ALT = -1 

DO 1530 1=1, TOUR 

IF!TCOST!I ) .GT.BEST!N,2) ) GO TO 1530 
ALT = ALT + 1 
1530 CONTINUE 

WRITE(6,I5) ALT 
WRITE<6,48) ITER 

1575 IF! ITER. LE . ITS) WRITE!6,55) TIMEX 
IF(ITER.GT.ITS) WRITE!6,56) TIMEX 

: WRITE OUT ITERATION INFORMATION 



WRITE(6,4) 

WRITE!6,5) 

IF!ITER.GT. ITS) ITER 
WRITE!6,31) 



ITER - 1 



DO 1600 1 = 1 , ITER 

16 00 WRITE!b,24) (YTAB! I , J ) , J = 1 , 7 ) , ! YBT A8 ( I , J ) , J= 1 , 2 ) 
RETURN 
END 
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SUBROUTINE MINELM 



IMPLICIT INTEGER’t^aiA-Z) 

REAL'^4 TIMEX 

INTEGER M(20,20, 19) ,THETA,MAXELtMINEL(20) ,MIN(20) T 
1TITLE( 17), ZG,WY 

DIMENSION A( 20,20, 19) , COST (20) ,LEGC0M(20) ,TC0ST(30) 
1FM(20) ,T0( 20) ,X(20) , ARCCOM ( 20 , 20 ) , BE ST (20,4) 
2YTAB(2500,7),YBTAB(2500,2) , I K ( 20 ) , JK ( 20 ) 

COMMCN TIMEX,M,THETA.MAXELtMINEL,MIN,TITLE,ZO,V;Y,A, 
1C0ST,LEGC0M,PM,T0,X, ARCCOM, BEST, YTAB,YBT A B, IK, JK, TOUR 
2AA,N,L , ALIKE, LEGREQ, STEP, I T ER , DEL , WX , MAX L EG , L EG , W YB AR 
310, JO,KG,ITS, TCOST, INDEX 



DO 630 K=1,L 
KEY = 0 
MINEL(K) = ~1 
IK(K) = 0 
JK(K) = 0 

IF (K.eQ.LEGCOM(K) ) GO TO 630 
602 MIN(K) = 32C00 
DO 629 1=1, N 
DO 629 J=1,N 
IF( I ,E0. J) GO TO 621 

IF ( K.EQ.l. 4ND.LEGC0M( 2) .EQ,0) GO TO 605 
IF (K,EQ.l ) GO TO 606 

IF (K,EQ.L.AN0.LEGCCM(L-1).EQ.0) GO TO 605 
IF (K.EQ.L) GO TO 610 

IF( K.EQ. LEGC0M(K+1 )-!• AND,K.NE.LEGC0M(K-1) +1 )GOTO 606 
I F( K. EQ. LEGCOM( K“1 ) +1 . AND. K. NE . LEGCOM ( K+ 1 ) -1 ) GOTO 610 
IF( K.E0.LEGCQM(K+1 )-l. AND. K. EQ . LEGCOM ( K- 1 ) +1 )GOTO 609 

605 IF (ARCCOM( I,J).LT.100) GO TO 621 
IF (KEY.EQ.O) GO TO 620 

IF (KEY.EQ.l) GO TO 622 

606 IF ( J.NE.FM(K+1 ) ) GO TO 621 

IF (ARCCOM( I,DEL).LT.I00) GO TO 621 
IF (KEY.EQ.O GO TO 620 
IF (KEY.EQ.l) GO TO 622 

609 IF( I .EQ.T0(K-1 ) .AND. J. EQ.FM( K4-1 ) . AND . KE Y . E 0. 0 ) G0T0620 
I F( I .EQ.TO( K-1 ) . AND. J. EQ.FM( K+1 ) .AND.KEY.EQ. 1 )GGT0622 
GO TO 621 

610 IF ( I.NE.TO(K-l) ) GO TO 621 

IF ( ARCC0M(DEL,J).LT.100) GO TO 621 
IF (KEY.EQ.l) GO TO 622 

620 MINEL(K) = M I NO ( MI N ( K ) , M ( I , J , K ) ) 

IF ( MINEL( K) .GE.MIN(K) ) GO TO 621 
IK(K) = I 

JK(K) = J 
MIN(K) = MINEL(K) 

621 IF (KEY.EQ.l) GO TO 629 
IF ( l.NE.N) GO TO 629 
IF ( J.NE.N) GO TO 629 

REDUCE MATRICES 



KEY = 1 
GO TO 602 

622 IF( I .EQ. J) GO TO 629 

M(I,J,K) = M(I,J,K) - MINEL(K) 

629 CONTINUE 

630 CONTINUE 



RETURN 

END 



SUBROUTINE CHECK 



C THIS SUBROUTINE DETERMINES IF A COMPLETE ROUTE CAN BE 
C ENUMERATED AFTER ONLY EVERY OTHER LEG HAS BEEN 
C DETERMINED, STARTING WITH LEG ONE 



IMPLICIT INTEGER*2(A-Z) 

REAL=^^4 TIMEX 

INTEGER^A M( 20,20, 19) , THETA , MAXEL , MI NEL ( 20 ) ,MIN(20) , 
1TITLE( 17) , ZG,WY 

DI MENS ION A (20,2 0, 19) , COST (20) , LEGCOM( 20) ,TC0ST(30) 
1FM(20) ,T0( 20) ,X(20) , ARCCOM ( 20 , 20 ) , BEST (20,4) 
2YTAB(2500, 7 ) , YBTAB ( 2500,2) ,IK(20),JK(20) 

COMMON TIMEX, M,THET A, MAXEL ,MI NEL, MIN, TITLE, ZO,WY, A, 

1 COST, LEGCOM,FM, TO, X, ARCCOM, BEST, YTAB, YBTAB, IK, JK, TOUR, 
2AA,N,L , ALIKE, LEGREQ, STEP, I TER, DEL, WX, MAXLEG, LEG,WYBAR, 
310, JO, KQ, ITS, TCO ST, INDEX 



DO 200 K=1,L,2 

IF (LEGCOM(K) .NE.O) GO TO 200 
RETURN 1 
200 CONTINUE 

LAST = N-4 
DO 300 K=2,LAST,2 
LEGCGM(K) = K 
FM(K) = TO(K-l) 

TO(K) = FM(K+1) 

300 COST(K) = A(FM(K) ,TO(K) ,K) 
RETURN 
END 
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SUBRCUTINE TIMEIT 



C NN=0 STARTS CLCCK: NN=-1 STOPS CLOCK 

IT = NN+2 
GO TO (20, 1C) t IT 
10 CALL TIMCN(f^M) 

TIMER = MM 
RETURN 

20 CALL T IMOEF(MM) 

TIME = MM 

TIME=(TIMEM-TIME) * 26.0 

RETURN 

END 



C ASSE 


NBLY LANGU/^GE LISTING TO CONDUCT 


TIMING ROUTINE 


TIMEALL 


CSECT 

ENTRY 


TIMON, TIMOFF 






TIMCN 


SAVE 


(14,12) ENTRY 


VIA 


-CALL TIMON(N)- 




USING 


T IMON, 12 








LR 


12,15 








ST 


12, TEMPI 








LA 


12, SAVEl 








L 


2, 0(1,0) 








L 


3,T0TIME 








ST 


3 ,CLOCKR 








ST 


3, 0(2,0) 








STIMER 


TASK,TUINTVL=CLOCKR 








L 


13, TEMPI 






EXIT 


RETURN 


( 14,12),T,RC=0 






TIMOFF 


SAVE 


(1A,12) ENTER 


VIA 


- CALL TIMOFF(N) 




USING 


TIMOFF, 12 








LR 


12,15 








ST 


13, TEMPI 








LA 


13, SAVEl 








L 


2, 0(1,0) 








TTIMER 


CANCEL 








ST 


0,0(2, 0) 








L 


13, TEMPI 








RETURN 


( 14,12),T,RC=0 








CNOP 


0,4 






TGTIME 


DC 


X'TFFFFFFF* 






CLCCKR 


DS 


F 






SAVEl 


DS 


IFF 






TEMPI 


DS 


F 








END 
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APPENDIX F 



MODIFICATIONS TO COMPUTER PROGRAM FOR 
DYNAMIC STORAGE ALLOCATION 

The following modifications to the basic program presented in 
Appendix E will provide dynamic storage allocation based on the number 
of nodes (N) and the maximum number of iterations (ITS) desired for each 
case in the computer data deck: 

1 . The Assembly Language listing on the following pages 
should be inserted at the very front of the computer source deck. 

2 . A new main program which is found after the Assembly 
Language listing replaces the original main program. The original main 
program becomes SUBROUTINE START and is listed here, after the new 
main program. 

3. All other subroutines remain the same as before with the 
exception of the variable type specification statements ^ DIMENSION 
statements, and COMMON statements. These 6 statements as found 
in the new SUBROUTINE START must be used in all the old FORTRAN 
subroutines except SUBROUTINE TIMEIT which does not change. 

4. The CALLS for the subroutines and the SUBROUTINE definition 
cards must be the same as before with the added arguments as found in 
the new SUBROUTINE START definition card. 

5. The JCL is included as a guide and is unique to the 

IBM 360 Model 67 Computer System installation at the Naval Postgraduate 
School. The only card which is required to be changed on various 
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runs in the EXEC card. It must contain a region for the GO step which 
is large enough to handle the case with the maximum number of iterations 
and specify time for the GO step large enough to accommodate the expected 
running time for all cases in the data deck. 

The makeup of the revised computer deck is on the following page. 
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INSERT DATA DECK 



C ASSE^^BLY LANGUAGE PROGRAM USED TQ OBTAIN DYNAMIC 
C ALLOCATION CF STORAGE SPACE BASED ON THE INPUT PARAMETERS 
C N AND ITS 





MACRO 

REGS 

LCLA 


6N 






LCLC 


CSYM 




.LOOP 


A MOP 






GSYM 


SETC 


* CN» 




R8SYM 


EQU 


CSYM 




6N 


SETA 


CN + 1 






AIF 


( CN LT 16 ) .LOOP 




MEND 






GETARY 


CSECT 

REGS 








STM 


RIA ,R12, 12IR13) SAVE REGS 




LR 


R12tR15 


GET BASE OUT OF R15 




USING 


GETARYtR12 






LA 


R11,SAVEAREA 






ST 


R13,4(R11) 


LINK SAVE AREAS 




ST 


R11,8(R13) 






LR 


R13,R11 . 






LR 


Rll.Rl SAVE ARGUMENT LIST LOC FROM GETMAIN 




USING 


ARGStRII 






L 


R2,ANUM 


GET NUM ARGUMENT 




L 


R3,0(R2) 


CHECK ITS VALIDITY 




LTR 


R3,R3 






BC 


13,NUMERR1 


NUM NEGATIVE OR ZERO 




LA 


R4,99 






CR 


R3,R4 






BH 


NUMERR2 


NUM GT 99 




LR 


RA,R3 


LENGTH OF SCRATCH = S’J^NUM + 8 




SLL 


RA,3 






LA 


R0,8(R4) 






GETMAIN R,LV=(0) 






LR 


R9,R1 






USING 


SCRATCH, R9 






ST 


R3,NUM START UNPACKING AND CHECKING CALL LIST 




L 


R5, ANEXT 






L 


R5,0(R5) 






CL 


R5, = X* OOOCOOOO 


* ADDRESS NEGATIVE CR GT 768K ? 




BH 


BADNEXT 






ST 


R5,NEXT 






SRL 


RA,1 ISSUE 


GETMAIN FOR NEXT ARGUMENT LIST 




LR 


RO,R4 






GETMAIN R,LV=(0) 






LR 


R1C,R1 






USING 


CALLLIST,R10 






LA 


R7,AL1 






LA 


R6,L1 






LA 


R4, ARRYl 




GETCORE 


L 


R5,0(R7) 


GET ADDRESS OF HUNK LENGTH 




LTR 


R5,R5 


LAST ONE ? 




BP 


CONTIN 


NOPE 




CH 


R3,=H‘1» DOES 


NUM AGREE ON LAST ONE? 




BE 


CCNTINl YES 






B 


NUMERR3 NO 




CONTIN 


CH 


R3,=H»1* DOES 


NUM SAY IT IS LAST AND IT ISNT? 




BE 


NUMERR4 




CONTINl 


L 


R5,0(R5) GET 


LNGTH 




CL 


R5,=X*00080000 


» LENGTH GT 512K OR NEGATIVE? 




BH 


BACLENG 






LR 


R0,R5 






GETMAIN R,LV=(0) 






ST 


R1 ,0(R4) 


PUT ADDRESS IN CALL LIST 




LA 


R4,4(R4) 


INCREMENT FOR NEXT LOOP 




ST 


R5,0(R6) PUT LTH IN SCRATCH FOR LATER FREEMAN 




ST 


R1,4(R6) 


PUT ADDRESS IN SCRATCH 




LA 


R6,8(R6) 


INCREMENT 




LA 


R7,4(R7) HUNK 


LENGTH AND SCRTCH AREA REGS 
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BCT 

SH 

MVI 

LR 

L 

BALR 

L 

SLL 
LA 



R3tGETC0RE 

RA,=H*4* 

0( P4) ,X*80* 
RItRIO 
R15,NEXT 
R1^,R15 
ROtNUM 
ROt2 

RItARRYI 



GET NEXT ARRAY 

PUT HEX 80 ON LAST ADDRESS 

PUT ADDRESS OR CALL LIST IN R1 

CALL NEXT ROUTINE 
GET RID OF CALL LIST 



FREECOR 



FREEMAIN R,LV=(0)tA=(1) 

LA R5tL1 LOOP TO FREE ARRAY CORE 

L R6,NUM 

L R0,0(R5) GET LENGTH OF FIRST ARRAY HUNK 

L R1,4(R5) GET ITS ADDRESS 

FREEMAIN R, LV= ( 0 ) , A= ( 1 ) 

LA R5,8(R*^) INCREMENT POINT TO HUNK LOC AND LGTH 
BCT R^tFREECOR 

L RO,NUM FREE SCRATCH AREA 

SLL R0t3 

LA ROtR(RO) NUM'^8 +8 IN RO 

LA RltSCRATCH 

FREEMAIN R,LV=(0) , A=( 1) 

L “ “ 



PROGRAM 





MVI 


0(P7) tXMO'BL 




MVC 


1(31,R7),0(R7 


RETURN 


L 


R13t4(R13) AL 




LM 


R14TR12t 12(R1 




SR 


R15tR15 




BR 


R14 


NUMERRl 


LA 


R6tMSGl 




B 


INSERT 


NUMERR2 


LA 


R6tMSG2 




B 


INSERT 


NUMERR3 


LA 


R6 ,MSG3 




B 


INSERT 


NUMERR4 


LA 


R6,MSG4 




B 


INSERT 


EACNEXT 


LA 


R6,MSG5 




B 


INSERT 


EADLENG 


LA 


R6,MSG6 




L 


R4tNUM 




LA 


R4,1(P4) 




SR 


R4tR3 numb 




CVD 


R4,MSG1 PUT I 




UNPK 


MSGl+8( 3) ,MSG 




01 


MSGU-IOtX ' FO* 




MVC 


MSG6+7( 2) , MSG 


INSERT 


L 


R7, AERRMSG 




MVC 


0(32tR7) ,0(R6 




B 


RETURN 




DS 


OF 


SAVEAREA 


DC 


18F*0* 




DS 


OD 


MSGl 


DC 


C* NUM OF ARR 


MSG2 


DC 


C* NUM OF ARR 


MSG3 


DC 


C* NUM OF ARR 


MSG4 


DC 


C* NUM OF ARR 


MSG5 


DC 


C* ADDR NEXT 


MSG6 


DC 


C* LENGTH 


ARCS 


DSECT 




AERRMSG 


DC 


A(C) ADD OF 


ANEXT 


DC 


A(C) ADD OF 


ANUM 


DC 


A(0) ADD OF 


ALl 


DC 


A(C) ADD OF 


ALNLM 


DC 


X*80« ,AL3(0) 


SCRATCH 


DSECT 




NEXT 


DC 


A(0) 


NUM 


DC 


A(0) 


LI 


DC 


F • C* 


A1 


DC 


A(C) 


LLAST 


DC 


F * C* 


ALAST 


DC 


A(C) 


CALLLIST 


DSECT 




ARRYl 


DC 


A(C) 


ARRYNUM 


DC 


A(C) 




END 





INSERT APPROPRIATE ERROR MSG 
AND RETURN DIRECTLY TO CALLI 
PROGRAM 



b(2) MSGl 

CHANGE BOTTOM ZONE TO FOX 
? PUT IT IN TEXT OF MSG6 



^YS - OR 0 
\YS > 99 

\YS > NUM GIVEN LNGTS 
< NUM GIVEN LNGTS 
ROUTINE - OR > 768K 
IS > OR > 512K 
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C NEW MAIN PROGRAM FOR DYNAMIC ALLOCATION OF STORAGE SPACE 



EXTERNAL START 
INTEGER^Z ALIKE, AA,NCASE 
INTEGER*^ EPROR(B) ,6LNK/4H / 
DIMENSION TITLE( 17 ) 

CCMMON/Z/ TITLE, N, ITS, AA, ALIKE 



1 FORMAT(IA) 

2 FORMAT! I^, 12,16, 17AA) 

25 FORMAT(8A4) 

READ(5,1) NCASE 

DO 10 AA=1, NCASE 

C READ INPUT PARAMETERS 

READ (5, 2) N, ALIKE, ITS, (TITLE! I ) ,1=1,17) 

C CALCULATE AMOUNT OF STORAGE REQUIRED FOR VARIOUS ARRAYS 

NNNM12 = N*N=<‘!N-1)'^2 

NNNM14=NN'NM12’<‘2 

N15 = 15 

N2 = N*'^2 

N4 = N^4 

NN2 = N=^N2 

N42 = N4=('2 

ITS72 = ITS^aA 

ITS22 = ITS’:«4 

C CALL ASSEMBLY LANGUAGE PROGRAM FOR OBTAINING STORAGE 

CALL GETARY! ERROR, START , N1 5 , NNNM 14 , N4, N4 , NNNM 12 , N 2 , N2 , 
1N2,N2,N2,NN2,N42 ,ITS72, 1TS22 ,N2,N2) 

IF(ERROR(l ) ,NE,BLNK) GO TO 20 
10 CONTINUE 

GO TO 30 

20 WRITE!6,25) ERROR 
30 STOP 
END 
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C THIS SUBROUTINE, START, IS THE MODIFIED OLD MAIN PROGRAM 



C' THE FIRST 11 CARDS OF THIS PROGRAM MUST BE USED IN ALL 
C OTHER FORTRAN SUBROUTINES EXCEPT SUBROUTINE TIMEIT, AND 
C THE CALLS AND CTHER SUBROUTINE CARDS MUST HAVE THE SAME 
C ARGUMENTS AS THOSE IN THE CALLS IN THIS SUBROUTINE, WITH 
C THE EXCEPTION CF THOSE WITH SPECIAL STATEMENT NUMBERS AS 
C ARGUMENTS IN THE OLD PROGRAM, AND THESE MUST APPEAR AT 
C THE HEAD OF THE ARGUMENT LIST. 



SUBROUTINE START ( M , M I NEL , M I N, A , COST , LEGCOM , FM , TO , X , 
lARCCCM, BEST ,YTAB,YBTAB , IK, JK) 



IMPLICIT INTEGER*2(A-Z) 

REAL*4 TIMEX 
INTEGER=i'4 N,ITS 

I NT EGER ={^4 THETA , MAXEL , T I TL E ( 17) ,ZO,WY 
I^TEGER=^'4 M(N,N,N) ,MINEL(N) ,MIN(N) 

DIMENSION A(N,N,N) ,COST(N) , L EGCOM ( N ) , TCOST ( 3 0 ) , FM ( N ) , 
1TC(N ) , X(N) , ARCCOM( N,N) ,BEST(N,4),IK(N),JK(N), 

2YTAB( ITS, 7) ,YBTAB( ITS, 4) 

COMMON TIMEX, THETA , MAXEL ,ZO,WY, TOUR, L ,LEGREO, STEP, DEL, 
IITER ,WX, MAXLEG,LEG, WYBAR, I 0 , JO , KO , TCOST , INDEX 
COMMON/Z/ TITLE, N, ITS, AA, ALIKE 



CALL INPUT (M,MINEL,MIN, A,COST,LEGCCM,FM,TG,X, 
1ARCC0M,BEST,YTAB,YBTAB, IK, JK) 

CALL ITERTE ( 82000 , M , M I NEL , M I N , A , COST , LEGCOM , FM, TO , X , 
lARCCOM, BEST,YTAB,YBTAB , IK, JK) 

CALL SOLN (M,MINEL, MIN, A, COST, LEGCOM, FM, TO, X, 
lARCCOM, BE ST, YTAB,YBTAB, IK, JK) 

2000 CONTINUE 
RETURN 
END 
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